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Abstract 


It is important to note that mixed systems of first and second-kind 
Volterra—Fredholm integral equations are ill-posed problems, so that solving 
discretized system of such problems has a lot of difficulties. We will apply 
the regularization method to convert this mixed system (ill-posed problem) to 
system of the second kind Volterra—Fredholm integral equations (well-posed 
problem). A numerical method based on Chebyshev wavelets is suggested 
for solving the obtained well-posed problem, and convergence analysis of the 
method is discussed. For showing efficiency of the method, some test prob- 
lems, for which the exact solution is known, are considered. 


Keywords: Mixed systems of first and second-kind Volterra—Fredholm in- 
tegral equations; Regularization method; Chebyshev wavelets; Convergence 
analysis. 


1 Introduction 


The Volterra-Fredholm integral equations [4,5, 7,13] arise from parabolic 
boundary value problems, from the mathematical modeling of the spatiotem- 
poral development of an epidemic, mathematical population dynamics, and 
from various physical and biological models. 
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In this paper, we consider mixed system of Volterra—-Fredholm integral 
equations (VFIEs) consisting of first and second-kind VFIEs as follows: 


t 
F(t,x) = AU(t, x) +f | K(t,7,2,5)U(n,s)dsdn, tel=(0,T], reEQ, 
0 JO 


(1) 
where F(t, 2) = [f(t,x), g(t, 2)]7, U(t,x) = [ur (t,x), ua(t, x)? and 
10 kii(t,7, 2,8) ky2(t,, x, 8) 
A= F i , K(t, 7,2, 8) = ae koo(t,, @, 8) 

Here, 2 denotes a (closed) bounded region in R? (d = 1,2,3) with the 
(piecewise) smooth boundary OQ. 

The reformulation of the initial-boundary-value problem for the linear 
heat equation in a two-dimensional spatial domain 2 with the boundary 0Q 
by single-layer techniques leads to a mixed system of Volterra—Fredholm in- 
tegral equations. Mixed systems of VFIEs are considered as the ill-posed 
problems. However, we will first apply the method of regularization that 
received a considerable amount of interest, especially in solving first kind 
integral equations. The method transforms mixed system to the system of 
second kind integral equations. The method of regularization was established 
independently by Phillips [12] and Tikhonov [14]. The method of regular- 
ization consists of replacing ill-posed problem by well-posed problem. Some 
numerical methods have been proposed for solving Volterra—Fredholm inte- 
gral equations of the second kind; see, for example, [2,6,8,10,11,15,16]. In 
this paper, the wavelet collocation method is developed for a mixed system 
of first and second kind Volterra—Fredholm integral equations. 

The paper is organized as follows. In Section 2, we consider some applica- 
tions of Volterra—Fredholm integral equation and in Section 3, we introduce 
regularization technique to transform mixed system (1) into the system of 
second kind integral equations. In Section 4, a numerical method based on 
Chebyshev wavelets is applied for solving the obtained well-posed problem. 
The convergence analysis is given in Section 5 and the numerical experiments 
are carried out in Section 6, which will be used to verify the theoretical re- 
sults. 


2 Some applications 


We can consider Volterra-Fredholm integral equation of the first kind as 
a spacial case of system (1)(let uw: = 0,ki2 = 0). The reformulation of 
the initial-boundary-value problem for the linear heat equation in a two- 
dimensional spatial domain 0 with boundary 02 by singlelayer techniques 
leads to a Volterra—Fredholm integral equation of the first kind in the follow- 
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ing form [3]: 


or(t,6) = | | K(t— s,X(0) — X(s))U(s, d)ddds, 


where X (0) is asmooth 1-periodic parametric representation of the boundary 
curve [ = OQ, and gr represents the function describing the given boundary 
condition on I x OQ. 


For the another example, you can consider the following integral equation 


1] 
f(a,t) = (A(t) — f(a) = i. [ F(t,7)k(, y)®(y, dydr, 


(x = E(r1,%2,£3),y = Y(Y15 Yo, Y3) (x,y) € Q,¢t € (0, 7], ) 
under the condition 


| ®(x,t)dx = P(t), 
fe) 


can be investigated from the contact problem of a rigid surface (G,v) hav- 
ing an elastic material occupying the domain 2 where f,(x) describing the 
surface of stamp. This stamp impressed into an elastic layer surface (plane) 
by a variable known force with respect to time P(t) whose eccentricity of 
application e(t), (t € [0,T]) that case rigid displacement A(t). Here G is the 
displacement magnitude and v is Poisson’s coefficient. 


3 Regularization technique 


The linear operator defined by second equation in system (1) has not gener- 
ally a continuous inverse, so that it is difficult to obtain a precise numerical 
solution by classical discretization methods. Thus regularization techniques 
could be used instead to transform integral equations such as second equa- 
tion in system (1) into second-kind integral equations. More precisely, we 
consider the following integral equations: 


t 
Fae: 7 | kia (t, 7,2, 8)ua(n, 8)dedn 
0 Q 
t 
+/ J batten s)ua.eln,s)dsin, 
0 fe 
= eeses | / kaa (t,9, 2, 8)u1 (9, 0)dedn 
0 Q 


t 
+/ i hoo (t, 1), X, S)U2,¢ (n, s)dsdn, 
0 7 


where € is a fixed positive number. 
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Under a hypothesis that will be clarified later, it can be proved that the 
solution uz-(#,t) of equation (2) converges (when ¢ — 0) to the solution 
Ug(x,t) of equation (1). Now, we consider the matrix form of equation (2) as 
the following system of integral equations: 


F(t,x) = EU,(t, x) +f [ Ktn2.5)Uclns)dsdn, (3) 
where E = F "| and U,(a,t) = [ui (a, t), us,.(a, t)]*. 


We need the following definition and lemma to regularization technique. 


Definition 1. A self-adjoint operator «: H + H, where H is a real Hilbert 
space, 1s called coercive if there exists a constant c > 0 such that 


(kz,2) > ellz||?, for all x € H. 


Lemma 1. Suppose that the integral operator of system (3) is continuous 
and coercive in a Hilbert space, where F', U, Uz are defined. Then 


e ||U-|| 1s bounded independently of €; 
e ||U- —U]|| tends to 0 when e > 0. 


Proof. From (3) we conclude the following: 


t t 
|Ell[ell = EF | [ KU-dsdn|| > —||F\| + | if fi KU-dsdn||, (4) 
0 


The coercivity property of the integral operator implies 


t 
| | i KU-dsdy]| > a|Uell, (5) 
0 Q 


where a is the coercivity constant. 
From (4) and (5) we deduce 


|E|\\|Uel] > al|Uell — ||F'll; (6) 
therefore it is obtained: 
(a — ||E\)) [Vell < Fl 


which proves the first part of the lemma. 
Now, for proving the second result, by using (1) and (3), we have 


EU. -U) =- i, - KGa. a, dt sy Ue Olde we 
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where U;(a,t) = [ui(2,t),0|". By rearranging (7), we can write 
t ~ 
| [ O22, 9)lUeln,s) — U(n,s)ldsdn = -EU.-0) 02, (8) 
0 JQ 
where, by using vectors operations EU — Ui — EU» is obtained, which U> — 


(0, u2]”. Taking the norm from the both side of (8) and using the coercivity 
property imply 


aU. — Ul| < || EU. — U) — eVg]| < ||B\\|U-e - Ul] + eal. 9) 
From ||U|| = ||u2||, we deduce that 
(a — ||E||)||Ue — UI] < elluall, (10) 


and now from ||£|| = 1+ ¢— 1 when e > 0 and a > 1 then it is concluded 
that ||U. — U|| > 0 when ¢ + 0. This completes the proof. 


The conditions of existence and uniqueness of solutions related to the 
VFIEs (3) can be investigated by considering the theorem about existence 
and uniqueness of solution of the second-kind Volterra—Fredholm integral 
equation in [3]. 


Theorem 1. Assume that 

1.FeEC(LlxQ) 

2.KeEC(D x?) where D = {(t,n),0<n<t< TH. 

Then the mixed system (3) possesses a unique solution U, € C(I x Q). 


4 Numerical treatment 


4.1 Chebyshev wavelets 


Dilations and translations of the “Mother function,” or “analyzing wavelet” 
®(x) define an orthogonal basis, our wavelet basis: 


®..1)(t) — 23 B(2-*¢ _ OF 


the variables s and / are integers that scale and dilate the mother function ® 
to generate wavelets, such as a Daubechies wavelet family. The scale index s 
indicates wavelet’s width, and the location index I gives its position. Notice 
that the mother functions are rescaled, or “dilated” by powers of two, and 
translated by integers. What makes wavelet bases especially interesting is the 
self-similarity caused by the scales and dilations. Once we know about the 
mother functions, we know everything about the basis. Chebyshev wavelets 
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P(nzm) = O(k, n,m, t) have four arguments, n = 1, 2,.. .,2*-1, where k can 
assume any positive integer and m is the degree of Chebyshev polynomials 
of the first kind. They are defined on the interval [0,1) as 


n 
Qk-1) 


(11) 


kes a Hi 
ete 


0, otherwise, 


where 
m= 0, 


a _ Te 


and m = 0,1,...,M—1, n = 1,2,...,2*-1. Also 7,,(t) are Chebyshev 
polynomials of degree m which are orthogonal with respect to the func- 


tion w(t) = —., on the interval [—1,1]. We consider the two-dimensional 
V1-@ 
Chebyshev wavelet $(n,m1,n2,m2)(£,y) as follows: 
27s Tn, (Qa — Ini +1) Tin, (2"y—2ng+1), Ba< ma 
a my ( z£— 2n4 + ) mo ( y — 2n2 + 5 ee 
phpat sy< pkg 1) 
0, otherwise, 
(12) 


where m, = 0,1,...,M, —1, m2 = 0,1,...,M2—1, m1; = 1,...,2"1—1, and 
mg =1,..., 227-1, 


A function u(x, y) € L?([0,1) x [0,1)) may be expanded as 


u(x, y) = > S- > S- Uny miname P(ni,mi,na,ma) (E> Y)- (13) 


ny=1m1=0 nga=1 m2=0 


If the infinite series in equation (13) is truncated, then 


2k1-1 y4,—12"2-1 Mo-1 


u(x, y) ~ es y S- S- Uni minome?(ni,mi nama) (Ls Y) = O(x)’UB(y). 


ny=1 m,=0 ne=1 m2=0 


(14) 
where ®(x) and ®(y) are (2*1~1)M, x 1 and (2*2—1) Mp x 1 matrices, respec- 
tively, such that 


(x) = [P(1,0) (Dlpetas 0(1,M,-1)(2), tees $(2'1-1,9) (2), oe) $(21-1,m,—1)(®)]" 


B(y) = [$(1,0)(y), - «+5 PC. Ma—1)(Y), «+ +5 Paka—2,0)(Y) «+s Pat2-1,My—1)(y)]- 


Also U is a (2*1~1)M, x (2*2-1) Ms matrix whose elements can be calculated 
as 


1 1 
benunee | | taaons) (2) aay U0, 8) Wor (2)0na (wave, 
0 0 
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where mz = 1,...,2"1-1, m, = 0,1,...,M@y —1, nz = 1,..., 2-1, me = 
0,1,...,M@M2-—1 and 


ky na-l ¢ n 
ee t-In+1), 33 <t< gh, 


otherwise 


Now, consider system (2) with J = Q = [0,1] and approximate the solu- 
tion of this system using (13) as 


ti (t, 2) = B(t)7U, G(a). (15) 


tiz,-(t, 2) = ®(t)? Us, ®(2). (16) 


Also, we approximate the kernels in system (2) respect to two variables 7 
and s as follows: 


Inserting (15), (16), (17), and the following collocation points 


a 


qe 12 12 MO, 
Moers ay eee 
J ieee 
“= oe) Pea a) 


into system (2), we have the following linear system of algebraic equations 
for the unknown coefficients U; and U2..: 


F(tis2;) = BL), 8( +f [ fia (tis, 23, 8)®(n)"Us®(s)dsdy 
+f i, hua(ti,n, 2), 8)®(n)? Ur,<(s)dsdn, 
g (ti, 2j) = €®(t;)" U2, -®(2;) +f [ea hex (tin, tj, 8)®(n)" U1 ®(s)dsdn 
“l [ faa(ti,n, aj, 8)(9)"U2,-®(s)dsdn, 


6=1,2,...,M\2"-1, 7 =1,2,..., Mo2'-1, 
(19) 
4.2 Normalization 


In the language of optimization theory, for the linear bounded operator K : 
X — Y and ka = y, determine x. that minimizes the Tikhonov functional 
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J-(x) = ||Kx — y||? + ellz||?, for alla e X. 


We can consider the following theorem from [9]. 


Theorem 2. Let K : X — Y be a linear bounded operator between Hilbert 
spaces and lete > 0. Then the Thikhonov functional Jz has a unique mini- 
mumaze € X. This minimum x, is the unique solution of the normal equation 


exe + K* Kaz = K*y. 
Here, K* :Y + X denotes the adjoint of K. 


By using Theorem 2, system (2) with J = Q = (0, 1] can be written in the 
normal form as follows: 


Dp 1 
(oui | | kis (p,7, @, 8% (n, 8)dsdn 
0 0 
p 1 
+f | ky (p, 7,49; 8)ua,e(n, s)dsdn, 
0 0 


1 1 
/ | koo(p,t, q, x) g(t, x)dadt 
p 0 
— EU2,¢(D, 


q) 

ab 1 t 1 

i | | koa(p, t, 4,2)ka1 (tn, @, 8) (n, 8\dsdndaedt 
0 0 0 

1 


1 t 1 
| | i koa(p, t, 4, 2)kaa(t, 2, 2, 8)t2,¢(n, s)dsdndedt, 
0 0 0 


(20) 
Now, we can consider the numerical method based on Chebyshev wavelets 
from the previous subsection for the approximate solution of system (20). 


+ 


+ 


I 
I 


5 Convergence analysis 


In this section, we investigate the convergence analysis of the proposed 
Chebyshev wavelet collocation method, using polynomial approximation the- 
ory. 


Lemma 2. Assume that u(a,y) € L?([0,1) x [0,1)) can be expanded in the 
form of series (13) and that ti(t,x) is the approximation of u(a,y) which is 
defined by (14). Then a(t,x) converges to u(t, x). 


Proof. We recall the series u(t,2) from (13) and the truncated series ti(t, x) 
from (14), respectively, as follows: 


oOoO.—lUCOUlmUmUmlUMmCUCUWD 


u(x, y) = S> > S- S- tas manama PCaxsinainajmaay es 9) 


ny=1m1=0 nga=1 m2=0 
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and 


2ki-1 4,—12"2-1 Mo-1 


= » y » » Unyminome2P(ni,m1,m2,ma) (ty 2), 


ny=1 m,=0 no=1 me=0 


where m, = 0,1,..., My, —1, m2 = 0,1,...,Me—-1, m = 1,2, 28-4, 
and ng = 1,...,2*2-1. The functions $(n, m1 ,n2,m2)(t,2) are the Chebyshev 
wavelet. Let L?({0,1) x [0,1)) be the Hilbert space and let 
ahi tne) at peur) 
Biri grits ria tian) (Es x)= = 27-2 e(2 2 t— m,)¢(2 a i mg), 
where $(n,,m,,n2,;m2)(t,£) form a basis of L?((0,1) x [0,1)) as in (12). 
From (13), we consider 


[oe Co 
u(t, x2) = S- S- hiniime? anda ee) Sy ae jy (t, x) 
m,=0 m2=0 1=0 7=0 


where Gj =< u(t,x), b¢,;)(t,2) > for ky = 1,kg = 1 and < .,. > rep- 
resents an inner product. Let us denote $(,;)(t,2) = @(t,x) and ay =< 


u(t, x), b(t, xv) >. 


Define the sequence of partial sums S,, 2 > m of {a;;o(ti,x;)}. Let 
Snri,m, and Sn5m, be arbitrary partial sums with nj > ng. We show that 
Sn m is a Cauchy sequence in a Hilbert space. 


Let Spr = 3 55 fe 2 9 ig O(ti, £;). From 


ny ne ni ne 
x), S- S- aij b(ti, rj) >= S> S- laa”, 
i=0 j=0 i=0 j=0 


2 ny n2 2 
we can show that |[Snjjn2 — Smymall” = Wiem, 41 j=me41 |Ois|” for ni > 
m1, ng > Moe. 
We can write 


ng 


2 
I SS ayotteznl 
t=m14+1 j=me24+1 


ng 


=< ¥ S a O( lay), ¥ S as P( tity) > 
=m, +1 7=mMmet+1 =m +1 7=mMe24+1 
ng 


> dD las, 


=m +1 7=me24+1 


for ny > m4,N2 > mg. Then it is concluded that 
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> 3 aig (ti, £5) IP = > 3 les), 


t=m1 +1 7=mMme24+1 =m +1 7=mMme24+1 


for ny > m4,N2 > M2. 
Now from Bessel’s inequality, we deduce that $7’ 


i=m4+1 ed |v |? 


convergent and therefore || )7y" 4.1 PBs mapa OG OCH By) ||? — 0, 

that is, || Tim .41 pphatare ai; O(t;,£;)|| + 0 and {Sm} is a Cauchy se- 
quence thus it converges to a real number like ‘s‘ 

In the continuance, we prove that u(t,x) = s. 


<s—ull,2), Oly, 24) > =< 8,0, a:) > — < ull,2), OG, a) > 
=< lim Snm, P(ti, Li) > Ai; 
noo 


= lim < Sang Oli Bi) > — Oe = By — ey = 0; 


then it is concluded that u(t,x) = s or 0/5 3742.5 aij O(ti, Zj) converges to 
u(t,xz), and by induction it is evident that this raeult is established for any 
integer numbers of kh; and ko, so for ky — oo and kg > oo. In the other hand, 
we can consider f(t, 2) as Sym, which means u(t, 7) converges to u(t, x) and 
this completes the proof. 


Theorem 3. Consider a function u(x, y) € L?({0,1] x [0,1]), with bounded 


4 
forth partial derivations, oe | < B. Then the wavelet coefficient, 


Uniminom, in (14) , decay as follows: 


TB 
Unimynem < * 
jtmiminamal S Sap ng) 8/2(ne — 1)(me = 1) 


(21) 


Proof. From (14) it follows that 


Uniminzm2 = [ [. (x iY Oni, mig) x) by (ne, my )Wny (x)Wn. (y)dydx 
=f Pn, mye Vite, (a (fw ula Y Ono, me) (Y )Wne (y)dy) dx. 


np /2*2-1 
= Pm, m)(£ Wr, (x Cf, 2h /Pa an, y) 


)/2k2-1 
Testor 2y — 2nz +1) 
Wn (2*2y — 2ne + 1) dy) dx. 
(22) 
If m1, m2 > 1, by substituting 2*2y — 2ng + 1 = cos@ in (22), it yields 
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1 : e cos@2ng—1, /2 
Qk2/2 [ (nq ,m1) (@)wny (x) (| u(x, —=ohe= =o . cos m28 do] 


Unjmingmg = 
v2 cos @ + 2n2 —1.,sinm20.|* 
~ ao (m1 ,my) (#) Wn (x) uc, ke ) as He 
2 . 6 — 1)/2k2 7 
tee = | ons ae / ) sini tngdiatneda dx 
2 4 
= So om itiod (yee 
Br oe ees L)Wry Dy 
sin(m2 — 16 sin(m2 + =) 7 
0) 


2 os 4 ee 
1 * 07u(x, (cos @ + 2ng — 1)/2*2) ] 
Rm (0) d6| de, 
sma | O2y (9) ‘liom 
(23) 


where 
sin(m2—1)@ — sin(m2 + 3 


Nm(8) = sin O( a ae 


Then, we obtain 
1 6? u(ax,(cos 0+2n2—1)/2*2) 
2 


1 Tv 
Yom = BT fig fy [Is 2 
Bama) ()WwE, () der] ng (8) dB. 


Now, similar to the discussion in (23), by substituting 2"*y—2n,; +1 = cosa 
in (24), it yields 
Unymyngmg 


1 


= 25(k1 +k2)/2m1mg(27) 
2n1 — 1)/2*1, (cos@ + 2ng — 1)/2*2 
cosa@ + 2n1 )/ > (cos a amg )/ ) hm, (a)hmg (8) da dd. 


™ f* du(( 
“ 2ra2 
o Jo 0? 207y 


Thus we have 


[Un ymingma| 


= 25(k1+k2)/2m mo2(2r) 
Qny — 1)/2*1 9 + 2nz — 1)/2*2 
cosa + 2n1 )/ , (cos ar ene )/ ) hm, (a)hmg (8) da dO 


aa O*u(( 
0 0 B = 02202y 
Vion (@)| a i [ma (8)| 49. 


~ 25(k1+k2)/2m1m(27) 
(25) 
However 
. "= i —1)0 i 1 
7 |Rms(0)| dd = | sin 9 (2mm i SeneT ,) \ao 
0 0 mg —1 ma+1 
™ sin @sin(m2 — 1)0 sin @ sin(m2 + 1)6 
< | | (26) 
0 m2 — 1 mez + 1 


2Mot 
~ (mz —1)’ 


138 S. Pishbin and J. Shokri 
and similarly, it is obtained 


2m, 


| [Rm (@)| da < (m2 —1)° (27) 


1 
Since n, < 2*1-1 and ng < 2*27!, by substituting (26) and (27) in (25), the 
desired result is obtained as follows: 


27B 
nimi N2M as . 
[Uniminamal S 55 +%)/2Gm2 — 1)(me — 1) 


Theorem 4. Let u(x,y) € L?((0, 1]x[0, 1]), with bounded forth partial deriva- 


tions, Ea | < B; then the error bound would be obtained as follows: 


Ok ,My,ko,M = Oe), (28) 
where 


Ok1,Mi,k2,M2 
2ki-1! 4,—-12"2-1 Mo-1 


-([ [ [uw uzy- > oY a tinymsnama Sina.) 2)P{nayma)(¥)] 


ny=1 m,=0 no=1 Mo 


L 


Wn, (2)Wn,(y) dx dy) . 


Proof. From the error statement, we have 


2 
Pk1,M1,k2,M2 
aki—1 My—129%2—1 Mo-1 


= [ [ [we,) ey be = b> Unymyngm2P(ny m1) (Z)P(na,ma) (Y y} 


ny=l m,=0 ng=1 M2=0 


Wn () Wns (y) da dy 


=f fe > 3 > Unymynam2P(n1, m4) (& )O(na, ma) (Y) 


ng=1 
aki—1 My—12"2—-1 Mo-1 


~ 


2 
Unimyngm2P(n4, m4) (& (no, ma) (Y ) Wn (2)Wny (y) dx dy 


iM 
iM 


ni=1 mi=0 ng=1 m2=0 


com co 


= [ P» 3 > tn seas (dni, m1) (® «)) ? (b(no,ma)(Y))” 


=2k1 my =M, ng=2k2 mg=M2 


Wn ()Wns (y) da dy 


co oo Co 


co 1 1 
S & & SY tis | / (CeCe 
10) (0) 


ny=2*1 m1=M1 ng=2k2 m2=M2 
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Wny (@)Wns (y) da dy 


ii x 2 
ny ‘1 k 
ay ane a aa jee [2 2 Tm, (2 ca a 2n1 ae | d 
= > ue. Tc 
1M1N2QM2 Ff ny-1 —_ (oki, _ 5) 
ny=2k1 ™1=M] ng=2k2 m2=M2 aki —1 vi (2-2 — 2m +1) 


ye oe ei 2 
pho ni [2 2 Tig (2"2y — 2n2 + | 
d 


x Yy 
— J1 — (2k2@ — Ing + 1)? 
a 2 
3 y y 3 gkitke Fi 7 [ie Qhe micas 2) d. 
a re ay ic 
ny=2ki ™1=M1 no=2k2 m2=Mg2 5hi= = Jl = (21a —2n, 4+ 1)? 
- 2 
Shoat Te (2%2y — 2ng + »] 
ak 1 — (222 — 2ne + 1)? 
Now, let Qhy —2n,+1=t; and Qhay — 2n2o +1 = to, then it is obtained 
ee) lee) lee) lee) 1 
3 2 7 ~ ~ ~ ~ 2 Tn (ti) (ti) 
SUGM Ok, My ,k2,M2 — oe s SS Unymineme 4 =e dty 
ny=2k1 m1 =My ng=2k2 M2=M2 
ae 
Ting (ta) dts 
a . 
-1VV1l-—t 
(29) 


For m, > 1, mg > 1, we have 


: Drm (t1) (4) 4 


TT 
1 Jin8" a fie 


then (29) simplifies as follows: 


rq co co co lo ) 
2 _ ~ ~ ~ ~ 2 
Ok,,M1,k2,M> — 4 ) ) ) ) Uniminems’ (30) 


ny=2k1 m=M, ng=2k2 M2=Mo2 


Therefore from (21) and (30), we can conclude the desired result as follows: 


1 
es = 3 3 3 Sy (nin2)>(m? — 1)2(m2 — 1)2° 


n= =21 miz= My, n2= =a m2= Mo 


(31) 
Also 
= 1 1 1 —— 1 
= a : 2 
a (n1)® (2'1)5 + (2' + 1)5 (2'1)5 = (1 4 oer)” (3 ) 


From (31) and (32), we conclude that 
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S35: 
Pte = o(2 re), 


Theorem 5. Let us consider M = max||K(t,n,2,s)||, for nt € I = 
(0, 1],s,¢ € Q = [0,1]. Assume that U.(t,x) is the exact solution of sys- 
tem (3) and that U.(t,x) denotes the Chebyshev wavelet approximation for 
the exact solution U which is given by (15) and (16). Then U-(t, x) converges 
to U-(t, x) and the following result can be obtained 

[Ue — Gell = O(2- B+). 


Proof. According to the proposed method in previous section, we consider 
(15), (16), and (17) and insert U-(t, x) and K(t, 7, x, s) as approximations of 
the exact solution U and kernel K into system (3) 


F(t,x) = EU-(t, 2) +f | K(t,, x, s)U-(n, s)dsdn. (33) 


Subtracting (3) from (33) and some manipulations, we get 


t 1 
zIllive — Gell < | | | Ce ee ACO LE | 
0 0 
t 1 2 7 
+I | | RY shan 8) Gay 0,0) dada: 
0 0 
From ||£|| > 1 it is obtained 
We — el] < |]K — Ky||Uel| +e — Gell. (34) 


Relation (34) together with Lemma 2, shows the convergence of the exact 
solution to the approximate solution. Considering Theorem 4, we have 


IK —K| =O Hered), (35) 
From (34) and (35), we conclude that 


We — Bel] = O(2- #4). 
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6 Numerical examples 


To demonstrate the efficiency and the practicability of the proposed method, 
we consider the following two examples. All results are computed by using a 
program written in the Mathematica®. 


Example 1. Consider the following TIAEs: 
t pl 25 
/ 7 100t(t + x)u(n, s)dsdn = P(B+41)(t+2), FE (0,1], (36) 
0 Jo 


where the exact solution is u(t,z) =t+ 2. 
Considering regularization techniques and (3) for (36), we have 


eu, (t, £) + [ [ 100t(t + x)u(n, s)dsdn = 2P(3 +4t)(t+2). (37) 


We assume that &-(t,x) is the approximation of the exact solution u-(t, x) 
which is defined by (15). For analyzing the behavior of the error representa- 
tions, we consider absolute error as 


Error = |ult;,%;) — at, 2;)|. 


The Chebyshev wavelet method described in Section 4 has been imple- 
mented for problem (37) with M, = Mz = 3,k, = kg = 2 and the error for 
different values of ¢ has been reported in Table 1, which confirms the theoret- 
ical results of Lemma 1. Figure 1 presents the plots of exact and approximate 
solution with ¢ = 0.000001, which are found to be in good agreement. 

From subsection 4.2, we consider the normal form of equation (36) as 
follows: 


1 pl pt pl 
EUe(D, g) + | | | 10000pt(p + g)(t + x)ue(n, s)dsdndadt 
p_Jo. Jo Jo 


2500 r 
a - 3 Pl + q)t*(3 + 4t)(t + x)dedt, pé€ (0, 1], 
p 0 


(38) 


and then solve this equation by the Chebyshev wavelet method with M, = 
Mz = 3 and ky = kg = 2. We report the error for different values of ¢ in Table 
2. From numerical results in Tables 1 and 2, we observe that the results in 
Table 2 are more accurate than the reported results in Table 1. 


Example 2. Consider the following: 


Au(t, x) +f | K(t,n, x, s)u(n, s)dsdn = f(t, x), t € [0,1], (39) 
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Table 1: Absolute errors of ue for different values of ¢ in Example 1 
Error Error Error 

(t;,«;) 2=0.01 2=0.001 == 0.000001 
(0.1, 0.1) 3.62 x 10-3 3.64 x 1074 3.65 x 10-7 
(0.2, 0.2) 1.62 x 10-4 1.45 x 1075 1.40 x 107-8 
(0.3, 0.3) 1.37 x 10-3 1.39 x 1074 1.39 x 10-7 
(0.4, 0.4) 1.32 x 1073 1.35 x 1077 1.37 x 10-7 
(0.5, 0.5) 9.09 x 1073 9.29 x 107 9.32 x 10-7 
(0.6, 0.6) 1.71 x 10-3 1.74 x 1074 1.72 x 10-7 
(0.7, 0.7) 1.69 x 10-3 1.73 x 1074 1.70 x 10-7 
(0.8, 0.8) 4.19 x 1075 9.27 x 10-7 1.87 x 10-9 
(0.9, 0.9) 7.42 x 10-3 8.20 x 1074 8.20 x 10-7 


Table 2: Absolute errors of ue for different values of ¢ for normal equation in Example 1 


Error Error Error 

(t;, 25) e = 0.01 € = 0.001 € = 0.000001 
(0.1, 0.1) 6.18 x 107 6.11 x 1O-* 1.95 x 1077 
(0.2, 0.2) 6.84 x 107° 6.26 x 10~° 522 10>? 
(0.3, 0.3) 3:79 10% 3.83 x 107-7 8.93 x 10-8 
(0.4, 0.4) DAs X10 2 R45 10°" 5.70 x 10-8 
(0.5, 0.5) 2.06:%.10-* 2.06 x 10~° 6.66 x 10-7 
(0.6, 0.6) 2.14 x 107° 214 xe 10% 2.42 x 10-7 
(0.7, 0.7) TIS LO Rie S105 * 1.67 x 1077 
(0.8, 0.8) 22510" 2.25 106% 2ALXIOT? 
(0.9, 0.9) Talo SLO? 7.18 x 10~® saa xA0TE 
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where 
{10 _ [P(s? +041) (s+) 
AS E AE AG MBE) = | 69 4 4 i? 8008 4 2008 | ? 
T od 
u(t, 2) = (u(é,2),ua(t,2)) , fe =(AGe. Ge) 4 


and f; and f> are such that the exact solution is 
wu(tjz)=?+2?+1, ue(t,z)=t+2. 


We apply the regularization method to convert the mixed system (39) to 
the system of the second kind integral equations. Then, the resulting second 
kind integral equation will be solved by the proposed numerical scheme in 
section 4 with M,; = Mz = 3 and kj = kp = 2. Let (tj, t%2,-) be the 
approximation of the exact solution (uw, u2,-) which are defined by (15) and 
(16). Numerical errors with several values of ¢ are displayed in Tables 3 and 
4. 


Figure 1: The plots of exact solution u and approximate solution of u with ¢ = 0.000001 
in Example 1. 


Similar to Example 2, we transform system (39) into the normal form of 
system (20) and then solve the normal system by Chebyshev wavelet method 
with M, = Mz = 3 and ky = kg = 2. The error for different values of ¢€ is 
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Table 3: Absolute errors of uy for different values of € in Example 2 
Error Error Error 
(i, y;) 2=0.01 2=0.001 == 0.000001 
(0.1, 0.1) 1.00 x 107-5 1.00 x 1076 1.00 x 1079 
(0.2, 0.2) 2.00 x 107° 2.00 x 10~® 2.00 x 10-9 
(0.3, 0.3) 3.03 x 107° 3.03 x 1076 3.03 x 1079 
(0.4, 0.4) 4.11 x 1075 4.11 x 107 4.11 x 1079 
(0.5, 0.5) 5.44 x 1075 5.44 x 1078 5.44 x 1079 
(0.6, 0.6) 6.58 x 107° 6.58 x 1076 6.58 x 1079 
(0.7, 0.7) 8.25 x 10-5 8.25 x 10-6 8.25 x 10-9 
(0.8, 0.8) 1.05 x 10-4 1.05 x 10-5 1.05 x 10-8 
(0.9, 0.9) 1.34 x 10-4 1.34 x 1075 1.34 x 1078 


Table 4: Absolute errors of u2,< for different values of ¢ in Example 2 


Error Error Error 


(xi, 4;) €=0.01 2 =0.001 = = 0.000001 


1.07 x 1073 1.83 x 1074 8.51 x 10-6 
1.58 x 1074 2.84 x 1075 1.86 x 107-5 
0.3, 0.3 2.39 x 10-4 2.36 x 1075 1.78 x 1075 
0.4, 0.4 4.63 x 1074 2.48 x 1075 2.37 x 10-6 


(0.1, 0.1) 
( ) 
( ) 
( ) 
(0.5, 0.5) 3.84 x 10-9 8.20 x 10~4 4.82 x 1075 
( ) 
( ) 
( ) 
( ) 


0.2, 0.2 


0.6, 0.6 1.66 x 1073 1.12 x 10~* 1.06 x 1074 
0.7, 0.7 122610 1.80 x 10~* 1.77 x 10-7 
0.8, 0.8 2.85 x 1073 2.81 x 1074 1.80 x 1074 
0.9, 0.9 6.69 x 10-3 3.88 x 107“ box 107? 
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Figure 2: The plots of exact solution ui and approximate solution of ui with ¢ = 0.000001 


in Example 2. 


Table 5: Absolute errors of uw for different values of ¢ for normal equation in Example 2 


Error Error Error 

(xi, yy) e=0.01 € = 0.001 € = 0.000001 
(0.1, 0.1) 8.54 x 1078 8.54 x 1079 8.54 x 10712 
(0.2, 0.2) 2.56 x 1078 256% 10-" 2.56 x 10-" 
(0.3, 0.3) 4.19 x 10-8 4.19 x 10-° 4.19 x 10-12 
(0.4, 0.4) 1.10 x 10-7 1.10 x 10-* 1.10 x 107" 
(0.5, 0.5) 5.46 x 1077 5.46 x 10% 5.46 x 10-14 
(0.6, 0.6) 1.73% 10° 1.73 x 10-7 Lyx 10- 
(0.7, 0.7) 2.54 x 10-5 2.54 x 107-7 2.54 x 10-1° 
(0.8, 0.8) 2.86 x 10~® 2.86 x 10-7 2.86 x 101° 
(0.9, 0.9) 2.58 x 10~® 2.58 x 10-7 2.58 x 10-19 


reported in Tables 5 and 6. Comparing displayed errors in the Tables 3, 4 
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Figure 3: The plots of exact solution uz and approximate solution of uz with ¢ = 0.000001 
in Example 2. 


and 5, 6, we observe that the results in Tables 5 and 6 are more accurate 
than the reported results in Tables 3 and 4. This is predicted by the theory, 
in particular by Theorem 2. Figures 2 and 3 present the plots of exact and 
approximate solution with ¢ = 0.000001. 


6.1 Haar wavelet method 


In literature, several wavelets with different properties have been derived 
and depending upon the applications, different wavelet families are used. 
The Cheybyshev wavelets method is found to be simple, efficient, accurate, 
and computationally attractive for solving linear and non-linear problems. 
The properties of Chebyshev wavelets make the wavelet coefficient matrices 
sparse which eventually leads to the sparsity of the coefficients matrix of the 
obtained system. The Haar wavelet is also the simplest possible wavelet. The 
technical disadvantage of the Haar wavelet is that it is not continuous, and 
therefore not differentiable. This property can, however, be an advantage for 
the analysis of signals with sudden transitions, such as monitoring of tool 
failure in machines. For comparison, we consider the Haar wavelet method 
to solve Example 1. The Haar wavelet family is 
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Table 6: Absolute errors of u2,¢ for different values of ¢ for normal equation in Example 


2 
Error Error Error 
(xi, Y;) é€=0.01 € = 0.001 € = 0.000001 
(0.1, 0.1) 1.25 x 1077 1.00 x 10-8 219% 10-7 
(0.2, 0.2) 9.92 x 10-7 9.83 x 1078 8.59 x 1078 
(0.3, 0.3) 9.58 x 1078 1.05 x 10-8 9.59 x 1078 
(0.4, 0.4) 3.64 x 1076 3.64 x 107-7 5.67 x 1079 
(0.5, 0.5) 5.75 x 1076 5.67 x 107-7 2.19% 107" 
(0.6, 0.6) 4.75 x 1076 A732 10-* 4.37 x 107-7 
(0.7, 0.7) 1.18 x 10-6 1,17 3 10-* 1,52 % 10-7 
(0.8, 0.8) 5.24 x 1076 5.24 x 107-7 7.34 x 1078 
(0.9, 0.9) 1.48 x 107° 1.48 x 10-6 1.39 x 10-7 
1, te [Fis Ta], 
h,(t) = —1, LE [72,73], 
0, elsewhere 
where 
k k+4 k+1 
fe 12S ? Cr 
m m m 
The integer m = 24,7 =0,1,...,p and k =0,1,...,m—1 are the level of 


the wavelet and translation parameter, respectively. The index 7 is calculated 
from the formula i = m+k-+1 and the maximal value is i = 2M(M = 2?). 
The index i = 1 corresponds to the scaling function of the Haar wavelet 
hi(t) = 1. We can consider an approximate solution of equation (37) as 


2M, 2M2 


ai.(t, x) = S- S- aighy(t)h; (x), 


i=1 j=l 


where the coefficients a;; can be obtained by inserting approximate solution 
a.(t,2) and the following collocation points in equation (37) 
i=) r—( 


i 
ty = (49... 3M, a=" eH 19 
l rh ee 14) ’ 1; x 2M> r 2 
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By assuming MM, = M2 = 4, the absolute error for different values of ¢ 
have been reported in Table 7. From Tables 1 and 7, we observe that the 
results obtained by the Chebyshev wavelet method are more accurate than 
the Haar wavelet method in this case. 


Table 7: Absolute errors of ue for different values of ¢ by Haar wavelet method in Example 
1 


Error Error Error 

(t;, 25) ée= 0.01 e = 0.001 € = 0.000001 
(0.1, 0.1) 8.99 x 107? 8.07 x 10-2 7.99 x 10-2 
(0.2, 0.2) 8.83 x 107? 211. 10-7 LAL 10-7 
(0.3, 0.3) AOA x 10-4 172 10-7 1.29 «10-7 
(0.4, 0.4) 9.33 x 10-7 7.99 x 10~? 7.82% 10-7 
(0.5, 0.5) 4.99 x 1071 5.00 x 107+ 4.99 x 10-4 
(0.6, 0.6) 9.58 x 10-2 8.25 x 10-7 8.08 x 107? 
(0.7, 0.7) 6.79 x 10-? 1.98. x 10-7 115% 10" 
(0.8, 0.8) 5.09 x 10-2 178 10-4 iW a a 
(0.9, 0.9) 9.26 x 10-2 8.01 x 107? T.84 «10-7 


7 Conclusion and future work 


This work has been concerned with the regularization method to convert the 
mixed systems of the first and second-kind Volterra—Fredholm integral equa- 
tions to the system of the second-kind Volterra—Fredholm integral equations. 
We presented the numerical method based on Chebyshev wavelets for solving 
the obtained second-kind problem. Convergence of the method was proved. 
These results were confirmed by some numerical examples. 


In the present work, we considered the mixed system (1) in special case 
with A = F 0| . The mixed system (1) in general form with A(t) will be 


investigated as our future work. 
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